Molecular heat pump for rotational states 
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In this work we investigate the theory for three different uni-directional population transfer 
schemes in trapped multilevel systems which can be utilized to cool molecular ions. The approach 
we use exploits the laser-induced coupling between the internal and motional degrees of freedom so 
that the internal state of a molecule can be mapped onto the motion of that molecule in an external 
trapping potential. By sympathetically cooling the translational motion back into its ground state 
the mapping process can be employed as part of a cooling scheme for molecular rotational levels. 
This step is achieved through a common mode involving a laser-cooled atom trapped alongside 
the molecule. For the coherent mapping we will focus on adiabatic passage techniques which may 
be expected to provide robust and efficient population transfers. By applying far-detuned chirped 
adiabatic rapid passage pulses we are able to achieve an efficiency of better than 98% for realistic 
parameters and including spontaneous emission. Even though our main focus is on cooling molecu- 
lar states, the analysis of the different adiabatic methods has general features which can be applied 
to atomic systems. 

PACS numbers: 33.80.Be, 37.10.Mn 
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I. INTRODUCTION 



The emerging field of cold molecules is a very vibrant 
topic in physics and physical chemistry. The consider- 
able interest in this topic is related to the properties 
of cold molecules and their many potential applications. 
Cold molecules have been identified as attractive systems 
for ultrahigh-resolution spectroscopy [l|, , quantum in- 
formation processing [3j], for developing new time stan- 
dards and testing fundamental physical theories such as 
the time variation of physical constants [J, [5l , the ex- 
istence of a dipole moment of the electron [fj, and for 
the measurement of parity violation Q- In chemistry, 
cold molecules are essential tools to explore quantum- 
mechanical effects in chemical reactions. In contrast to 
atoms, molecules have a very complicated level structure 
that consists of vibrational and rotational states as well 
as electronic levels. This abundance of states is the main 
obstacle for direct laser cooling of molecules. Usually, 
molecules do not provide the closed transitions required 
for cooling and non-destructive state-selective detection. 
This makes it impossible to perform direct spectroscopic 
measurements on single molecules — a standard technique 
in atomic physics. Additional complications result from 
the small energy differences between the rotational lev- 
els, leading to a thermal distribution of the population 
over the molecules' ro-vibrational states. 

Despite important achievements pj-IIol|. the control of 
molecular states never caught up with that of atomic 
systems. However, there has been remarkable progress 
in the synthesis of ultra-cold alkali dimers from sam- 
ples of ultra-cold atoms; see e.g., Refs. [IlT - fl3j . Fur- 
thermore, methods which enable the preparation of more 
diverse (e.g. polyatomic) cold molecular species in their 
vibrational ground-states have been successfully demon- 
strated. These methods include: supersonic beam ex- 



pansion followed by Stark deceleration [9J, optical Stark 
deceleration [lj| , electrostatic velocity selection [l4[ , col- 
lisional cooling in crossed molecular beams [l5[ and buffer 
gas cooling [8[. The vibrationally cold (but rotationally 
hot) states that result will be taken to be the starting 
point for the schemes described in this paper. The ex- 
perimental advances which have enabled the production 
of these cold molecular states have inspired theoretical 
investigations of the cooling of molecules by laser pulses 
[I6H2OJ or even by coupling molecules to an optical cav- 
ity [21|. Bartana et al. [18| used the electronic excited- 
state as a heat reservoir in order to cool the vibrational 
states of the electronic ground-state by means of short, 
shaped, laser pulses. In later work [19[ they employed 
state selective optical pumping, hiding the target state 
in a dynamically trapped state. Through this, Bartana 
et al. achieved a vibrational ground state population of 
97% after only 25 vibrati on p eriods. A related scheme, 
investigated by Schirmer [16| . increased the vibrational 
ground-state population to a similar level. These efforts 
to cool the internal degrees of freedom focus on the widely 
spaced vibrational states of molecules. In Ref. (2(| Bar- 
tana et al. investigated the possibility of cooling the rota- 
tional degrees of freedom in a simplified model employing 
the same techniques. However, even though the results 
of their calculations are very promising, the model has 
limited application as it neglects the vibrational degree 
of freedom of the molecule. 

In contrast to neutral molecules, ionized molecules can 
be sympathetically cooled by trapping them alongside 
atomic ions in a Paul trap. Under these conditions, state 
sensitive ultra-cold chemical reactions have been mea- 
sured [H, [23| and high resolution spectroscopy has been 
demonstrated on small ensembles (24|. Despite these 
achievements, the rotational degrees of freedom could not 
be controlled, but they led to new laser cooling schemes 
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for the internal states of molecules [ID, [26| exploiting the 
unique properties of these systems. 

In the present paper, we show that by means of purely 
coherent manipulations of internal states (i.e. rotational 
states), and by using sympathetic cooling of motional 
states, single molecular ions can be cooled close to their 
motional and rotational ground-state. The internal vi- 
bration of the molecule needs to be initially cold, in or- 
der for the final state to be cold in all its degrees of 
freedom. This initial state can be achieved using the 
existing methods described above. Then the cooling of 
the internal molecular state is achieved in three steps: 
first a laser cooled atomic ion is trapped alongside the 
molecule and a common mode of vibration is used to 
prepare the molecule in its motional ground-state. Next, 
using adiabatic passage methods |27H36| . the thermal 
internal state of the molecule can be mapped onto the 
molecule's motion in an external trapping field with high 
fidelity. During this mapping process the internal (rota- 
tional) states are transferred to the ground-state whereas 
the molecule's motion is excited. Finally, the molecule's 
motion is sympathetically cooled back into its ground- 
state using the common mode with a cooled atomic ion. 
By doing this without exciting other degrees of freedom, 
a molecule which is vibrationally, rotationally, and trans- 
lationally cold can be obtained. The overall process is a 
kind of molecular "heat pump" : the heat energy in the 
rotational degree of freedom is transferred by the coher- 
ent processes to the motional degree of freedom. This 
heat energy is in turn transferred to the environment by 
means of conventional cooling techniques, e.g. side-band 
cooling: this ensures that the process is uni-directional. 

Vogelius et al. [25[ investigated the cooling of molecular 
ions by coupling a single rotational state to the motion 
of the ion. By means of black-body radiation the popula- 
tion is pumped into the rotational ground-state resulting 
in a ground-state population of about 80% after a cooling 
time of the order of minutes. In this paper we focus on the 
cooling process by employing techniques from coherent 
control providing much faster and more efficient cooling. 
The result is a robust and highly efficient cooling process 
which enables the deterministic manipulation of the in- 
ternal states of molecules. We will examine adiabatic pas- 
sage schemes which exhibit high fidelity in conjunction 
with relaxed requirements on the experimental parame- 
ters compared to direct Raman transitions. With cooling 
times of the order of milliseconds, and final ground-state 
populations of more than 92%, the proposed scheme pro- 
vides a fast and efficient method for preparing molecules 
in their ro- vibrational ground-state. 

Even though we focus here on the cooling of molecular 
rotational levels, the technique is also applicable to other 
multi-level systems. So the technique can also be em- 
ployed in atoms with complicated level schemes, or more 
general ro-vibrational states of molecules. 

The paper is organized as follows: in Section [Til WG 
present the model used for our calculations. To find 
the best adiabatic passage process for our application we 



compare the results of numerical simulations for stimu- 
lated Raman adiabatic passage (STIRAP), Stark chirped 
Raman adiabatic passage (SCRAP) and chirped adia- 
batic rapid passage (CARP) in a A-type level system in 
Section HTT1 We also choose parameter ranges which are 
relevant for a possible experimental implementation. Sec- 
tion |IV] contains the results for an extended level scheme, 
and in Section [V] we present our conclusions. 



II. THE MODEL 

The model we employ in this paper is based on the 
states of a quantum mechanical rigid rotator which is 
a good approximation for the rotational states of small 
diatomic molecules. However, the techniques described 
here are applicable to most other level structures with the 
sole requirement that allowed Raman transitions between 
the states involved exist. In general, the method can be 
applied to ro-vibrational states of molecules as well as 
Zeeman and hyperfine levels. In order to simplify the 
discussion, and to avoid specializing to molecules with 
specific symmetries, we will not apply specific selection 
rules to the Raman transitions involved. Nevertheless, 
the results presented in this paper can be applied to a 
particular system by imposing the specific selection rules 
for that case with an appropriate relabelling of states. 
For example, we could utilize A J = 0, ±2 for linear 
molecules. 

We take the energy Ej of the rotational levels of the 
electronic ground state to be Ej = B ■ J( J + 1), with the 
rotational quantum number J and a rotational constant 
B. In order to limit the number of levels in our calcula- 
tions, only rotational levels up to a cut-off are considered, 
i.e. J < Jmax- For typical, light, diatomic molecules at 
room temperature only rotational states with J < 20 
are significantly populated. This decreases to below 10 
states for rotational temperatures lower than about 50K. 
This kind of temperature can be easily achieved by su- 
personic beam expansion [37| . In our scheme the levels 
J are coupled by laser pulses to an electronically excited 
state |e) (see Fig.[T]). From the excited state the molecule 
can spontaneously decay back into the electronic ground- 
state with a rate Tj. To represent the decay of the ex- 
cited state into levels outside the system, e.g. rotational 
states with J > J max or vibrationally excited states, an 
additional level \u) is also included in the model. The 
excited state can decay into this uncoupled state with a 
rate T u . 

The single molecule is trapped in a harmonic potential 
as provided by rf-traps for molecular ions. The quan- 
tized motion of the molecule in the trapping potential 
gives rise to an equally spaced ladder of motional states 
in addition to the internal states of the molecule. To take 
the molecule's motion into account, we use the standard 
notation with i representing the states J), |e) and 
| it). The quantum number describing the motional state 
is n. As discussed in the Introduction, in order to prepare 
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j = o 



FIG. 1: (Color online) Overview of the molecular levels and 
the couplings which are included in the model. Near-resonant 
laser pulses are driving the transitions between the excited 
levels \e,n) and the rotational states \ J,n). The excited state 
decays towards the states |J) at a rate Fj and towards the 
uncoupled state it) at a rate F u . The laser pulses have time- 
dependent Rabi frequencies Q a j p and are detuned from the 



relevant transition frequency by A' 



the molecule's motion in its ground-state, it is trapped 
alongside one or more atomic ions which can be directly 
laser cooled. The two types of ion form a crystal-like 
structure due to their mutual Coulomb repulsion, which 
enables sympathetic cooling of the molecule. In experi- 
ments with two types of atomic ion, ground-state popu- 
lations of better than 95% have been achieved [331. ^ n 
the molecular case, the minimal system would be one 
trapped molecular ion and one trapped atomic ion. In 
this case, we find that due to the frequency splitting of 
the motional modes, one mode (COM or stretch mode) 
can be singled out and used for the proposed scheme. 
The other mode only imposes an additional limit on the 
laser pulse length as we will discuss later in this paper. 

The laser-molecule interaction, in the rotating wave 
approximation (39| . is described by the following Hamil- 
tonian 1401: 



H(t) =hva^a + ]T HAP(t)\e)(e\ 
,/=i 

+ ^>[A^)-A^) + H|J-1)(J-1| (1) 
j=i 

+ [H k ;{t) + H k r {t) + H^{t)_ , 



where H*(t) represents the carrier resonance [4l[ for ei- 
ther the pump (k = p), or the Stokes-pulse (k = s), with 



J=l 



(2) 



O I + <7_ 



J=l 



The corresponding red sideband transitions in Eq. (fTJ) are 
the H^(t), which are given by 



H ?W = 2^ o [cr+a + a^a 1 j, 



J=i 



K(t)=J2 



(3) 



.7 = 1 



and for the first blue sideband transitions the H k {t) are 
given by 



H bW = }^ — o — V + +(7 - a )' 

.7=1 

j=i 



(4) 



In Eqs. (HE]); the secular frequency of the molecule in 
the external trapping potential is v and 77J are the cor- 



responding Lamb-Dicke parameters [41|. The detuning 
between the k-th laser frequency ui k j and the transition 
frequency w eJ of the \J) — > |e) transition is Aj(/j) = 
uj eJ — uo p ,{t) and A s T (t) = uj cJ _ 1 — uij(t) — v. The rais- 
ing operator for the internal states is ai_ = \e)(J\ and 
the lowering operator is a J _ = \J)(e\. The creation and 
annihilation operators for the motional number states 
\n) are a^, and a respectively. Note that wc work in an 
interaction representation with explicit time-dependence 
removed, and keep both the resonant and non-resonant 
couplings. The reason for this is that in pursuing the adi- 
abatic limit in Section Hill we will consider "strong" Rabi 
frequencies (ftj ~ v) which do not allow us to make a 
second RWA on the Hamiltonians for the sideband tran- 
sitions. 

We start our calculations with a molecule in an inter- 
nal thermal state such that it is already cooled in its in- 
ternal vibrational mode (e.g. by the methods mentioned 
in the Introduction) and such that the excited state |e) 
is not populated. We also assume that the molecule's 
vibrational motion in the trap has been cooled (e.g. by 
sympathetic sideband cooling [38[ ) so that only the man- 
ifold of rotational states | J, 0) are populated. The density 
matrix of this initial state is given by: 



fc=s,p 



Pinit = ~~ ^ (2 J 
J=0 



i) e -/^(J+i)| J)0)(0)J | ; (5 ) 
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where J m ax is the cut-off introduced for the numericai 
calculations. The normalization factor Z is given by Z = 
Ej=" (2 J + 1) e - 0- BJ ( J+1 ) with /3 = l/(k B T), and with 
T as the internal rotational temperature of the molecule. 

Starting from such an initial distribution, we will ap- 
ply coherent control techniques to transfer population 
between the different states. We aim to have a state 
mapping of the form 

£ Pj,o\J,0)(0, J\ -> £ fo,„|0,n)(n,0|, (6) 

,7=0 n=0 

where Pjo = Po.n = e~ /3BJI - J+1 ' ) /Z are the populations 
of the initial states J, n = 0} and the final target states 
| J = 0, n) respectively. A sequence of pulses will be used 
to map population in each J-state to a corresponding 
n-state with J = 0. 

Throughout this work, we assume that the system is 
initially prepared in a state given by Eq. ([5]), and de- 
rive the requirements for achieving the state mapping in 
Eq. ©, i.e. after completion of a number of the coherent 
pulse sequences. Thus we obtain a superposition of just 
the motional states which can then be cooled to the mo- 
tional ground-state \ J = 0, n = 0) by applying the sym- 
pathetic cooling [231 . [The state mapping @ can also 
be employed for non-destructive state detection: e.g. by 
measuring the initial thermal distribution ([5]). By cou- 
pling the electronic state of an atom trapped alongside 
the molecules to its motion, the mapped state can be 
read out.] For the coherent mapping the key idea is to 
use pairs of pulses, 5 (t) and ^ij(i), to induce population 
transfer between the states | J, n) and | J — 1, n + 1), see 
Fig. [TJ For this step it is important to have a resonance, 
so that additional states do not get strongly involved and 
disturb the mapping. Here, the resonance is arranged so 
that the quantity J + n is conserved at each step. (Al- 
though this is the simplest way to make the mapping, it 
would not be the only way, as we only require the transfer 
of population between unique pairs of states.) 

Repeating the J, n — > J — 1, n + 1 step J m ax times, 
where in each step J is decreased by one, the distribution 
of population can be moved to the J = motional states 
as in Eq. <j6j) . Since states with n are involved in 
the intermediate steps, it is clear that if we do not start 
in the motional ground-state the final state need not be 
entirely J = 0. However, when we start in the motional 
ground-state (i.e. n = 0), the population transfer is uni- 
directional. Hence, the population is transferred solely 
to the lower lying rotational states. 

In our analysis we numerically integrate the master 
equation for the density matrix p(t): 

^ = -mt),p{t)]+i(p{t))+u P {t)). (?) 

The Hamiltonian H(t) is given in Eq. ([TJ. The two Liou- 
ville terms C(p(t)) and C u (p(t)) describe the irreversible 
decay of the excited electronic state |e) to the rotational 



levels | J) of the electronic and vibrational ground-state 

,7=0 

(8) 

and the decay out of the system, i.e. to the uncoupled 
state | it), 

C u (p(t)) = -^-(p(t)alal+alalp(t)-2alp(t)aiy 

(9) 

In Eq. ^ the raising and lowering operators <j\ have the 
same form as o^, with the substitution J) — > \u). 

We will assume that the laser pulses are Gaussian, with 
a fixed width T and Rabi frequencies 

(10) 

where J = 1 , • • • , J max ■ The delay between the two 
pulses fij(i) and Q P j(t) that drive the transition | J) — > 
| J — 1) is 2t, whereas the delay between the J and J + 1 
pulse pair is f . The corresponding detunings will be ei- 
ther constants or have the form of time-dependent fre- 
quency chirps. We take all the Lamb-Dickc parameters 
to be the same, i.e. i] k j — r/, and for simplicity we assume 
that the decay rates are the same, i.e. r, 7 = r„ = T. In 
order to compare the three methods (STIRAP, SCRAP 
and CARP) in Section IIIII we characterize the popula- 
tion transfer efficiency by a parameter e representing the 
total population of the rotational ground-state after the 
transfer: 

e= A),n(* = oo), (11) 

n=0 

where Po,n(i = oo) is the population of the |J = 0, n) 
state after the transfer. For an efficient state mapping, 
the ground-state population (J = 0) will have increased. 
However, if the initial state also has some population 
in the |J = 0,n = 0) state, the efficiency e may also 
decrease due to laser- induced transfer out of | J = 0,n = 
0). Thus the definition (jll[) is not only a measure for 
the transfer efficiency, but also for the uni-directionality 
of the mapping process. For ideal state mapping the 
efficiency measure reaches the limit e = 1. In order to 
test the various passage methods in the next Section, the 
initial state is taken to be a mixture of 70% rotationally 
excited states (J > 0) and 30% ground-state (J = 0) 
population for each adiabatic method. This approach 
will test how uni-dircctional the scheme is. 



III. ADIABATIC PASSAGE METHODS 

To transfer the population from the state | J, n) to 
| J— l,n+l), various coherent processes can be employed. 
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FIG. 2: (Color online) The basic A-configuration used to in- 
vestigate the transfer efficiency of various adiabatic passage 
schemes in Section [IIII The pump and Stokes pulses, fif(t) 
and Oq {t) respectively, will induce a two-photon Raman tran- 
sition (which may be chirped) from state \n = 0, J = 1) to 
state | n = 1, J = 0). The best results are obtained if the 
initial population in \n = 0, J = 0) remains there. The off- 
resonant intermediate states |e, 0) decays towards the two ro- 
tational states and towards the uncoupled state |u) at a rate 
r n . The decay rate Fj represents decay from the levels e 
to the levels n, J. Other off-resonant states are included as 
shown. 



Here we focus on stimulated Raman adiabatic passage 
(STIRAP) dHHIH, Stark chir P ed Raman adiabatic 
passage (SCRAP or SIARP) 
batic rapid passage (CARP) 
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and chirped adia- 
which offer highly 
efficient population transfer in combination with robust- 
ness against variations of the pulse parameters. In order 
to simplify the numerical simulations we first investigate 
these processes in a system with just two rotational lev- 
els. That is, we examine in detail a single step in our 
multi-pulse coherent transfer scheme. In this case the 
most important states are the two lowest motional states 
for J = and the lowest (n = 0) motional state for J = 1 
(see Fig. However, to include the off-resonant effects 
all nine of the states shown in Fig. [5] are included in 
the numerical calculation. We choose parameter ranges 
which are relevant for an experimental implementation. 
The electronically excited state can decay to both rota- 
tional states as well as into the uncoupled state \u). 



A. Stimulated Raman adiabatic passage 



STIRAP is widely used in the optical control of 
molecules [HI, [H, H|[ where a Stokes pulse and a pump 
pulse are used in a "counter-intuitive" order to transfer 
the population between two states. The main require- 
ment for STIRAP is the two-photon resonance condition 
which corresponds to A^(i) = Af (t) for a transition on 




FIG. 3: (Color online) (a) The population transfer for STI- 
RAP with different delay times r and pulse widths T. The 
calculations are for the A-system in Fig. [2] The Rabi fre- 
quency is Qq = 1//20. (b) The efficiency as a function of 
Rabi frequency for T = 4500/^, r = S500/u. In (a) and 
(b) the pump and Stokes pulse detunings are A± — Af = 
and the other parameters are: Tj = r„ = 0.01^, r\ — 0.1, 
Po.o(-oo) = 0.3 and P li0 (-oo) = 0.7. 



a motional sideband (see Fig. [I}. In addition to this the 
pulse area must be large, floT 3> 1, and the delay be- 
tween the Stokes and pump pulse must be of the order 
of the pulse width, r rj T, to ensure the adiabatic evo- 
lution of the system. Another constraint arises from the 
necessity to address particular motional sidebands. The 
narrow splitting of the motional states imposes the use 
of laser pulses with narrow bandwidth vT 3> I [3(| . Fast 
pulses will inevitably result in the coupling of the tar- 
get state to other, close-lying states. This in turn will 
reduce the transfer efficiency. Furthermore, the resolved 
sideband condition fill ] requires that the Rabi frequencies 
are small compared to the trap frequency, i.e. Qq «C v. 
Thus, only slow and weak pulses can be used which, as 
we will see, substantially reduces the efficiency of STI- 
RAP. However, this can be overcome by employing other 
adiabatic passage schemes as described below. 



In Fig. 3(a) the transfer efficiency e is plotted for dif- 
ferent pulse widths T and delay times r. As mentioned 
above the initial mixed state is described by the popula- 
tions Po,o(~ oo ) = 0-3 and P\fi{— oo) = 0.7. The chosen 
decay rate of T = 0.01 v is in the typical range of values 
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for the decay rate of an electronically excited molecule, 
when compared to a typical trap frequency v (of the order 
of several MHz). The results in Fig. 3(a) show that the 



efficiency e is below 76%. The performance of STIRAP 
for short pulses is relatively poor due to the violation of 
the adiabaticity requirement (QqT 3> 1) and the limi- 
tation on the Rabi frequency fio posed by the resolved 
sideband condition (fio "C v). For long pulses the effi- 
ciency of STIRAP is compromised by the spontaneous 
decay r„ j, so the compensation of small Rabi frequen- 
cies by long laser pulses is not an option. Therefore, the 
efficiency of STIRAP in this parameter range is low and 
the population transfer is governed by optical pumping 
rather than coherent evolution. This is particularly vis- 
ible in the plateau region for t»T. In this regime the 
pulse delay is too large to sustain the adiabatic evolution 
of the system, which results in a net loss of the ground 
state population. Highly efficient, fast STIRAP between 
motional states requires large motional frequencies which 
are beyond current ion-trap technology. 

We can try to suppress the excited-state population by 
detuning the Raman transition from the excited state. 
However, it has been known for some time that detuning 
adversely affects the STIRAP process in the absence of 
decay j45(. With decay present one has to consider the 
balance of the adverse effect of detuning against a pos- 
sible reduction in spontaneous emission from the excited 
state of a model A system. Studies with such systems 
support the suggestion that the minimal losses (for mod- 
erate decay rates) are found by remaining on resonance 
[H |H Eg. Figure H shows how the efficiency of STI- 
RAP drops, for our model system, as we detune from 
resonance. The parameters are those for Fig. |3(b)| with 
the Rabi frequency ilo chosen to be at the peak of effi- 
ciency m Fig. [3(b)] We see that both with, and without, 
decay processes it is best to be resonant. In the case 
Tj = the resonance is much sharper, however. 

To exploit the adiabatic evolution of the system and 
to suppress the excited-state population by far detuning 
the Raman transition leads us to Stark chirped Raman 
adiabatic passage (SCRAP) and chirped adiabatic rapid 
passage (CARP), which will be discussed in Sections !!!! Bl 
and III! Cl In the limit of far detuning, the excited state 
population is strongly suppressed and the system's dy- 
namics are effectively that of a two-level system [48| . The 
effective Raman coupling between the states | J, n) and 
| J — 1, n + 1) is 



A5W ; 



(12) 



and the effective splitting of the coupled rotational levels 



A J (t) = 5 J (t) + S%t)-S p ,(t), 



(13) 



with the effective two-photon Raman detuning Sj(t) = 
A s ,(t) - A p j(t), and the two Stark shifts S s ,{t) = 
[7?0*(i)] 2 /4A*(i) and S p ,{t) = [n p (t)} 2 /4A p (t) induced 
by the Stokes and pump pulse respectively. 




FIG. 4: (Color online) Efficiency of the STIRAP process as 
a function of detuning (= Af). We show both the case 
with decay, Tj — T u — O.Olv (solid line), and without decay, 
Tj = r u = (dashed line). The other parameters are as in 
Fig.[3]with the best values taken for T = 4500/;/, r = 3500/is, 
and fio — v/20. 



Within this effective two-level model, adiabatic rapid 
passage techniques (ARP) [3(| can be applied. The main 
idea behind ARP is to drive the system through the reso- 
nance (Aj = 0) adiabatically, to achieve a complete pop- 
ulation transfer. The technique of Stark chirped Raman 
adiabatic passage (SCRAP) j3lTj33l | takes advantage of 
the Stark shifts whereas the chirped adiabatic rapid pas- 
sage (CARP) {27j - l3"o| uses overlapping laser pulses along 
with frequency chirps to transfer the population. We 
turn to these methods in the next sections. 



B. Stark chirped Raman adiabatic passage 

For the case of SCRAP (earlier known as self-induced 
adiabatic passage, or SIARP [3l|) the laser pulses are en- 
gineered so that the system undergoes an avoided level 
crossing (A,/ = 0) near a maximum of the effective cou- 
pling Qj(t) induced by the Stark shifts due to the delay 
of the pump and Stokes pulses. In order to obtain an 
efficient population transfer the system needs to evolve 
adiabatically in the crossing region [331 ]. For Gaussian 
pulses this leads to the condition: 



T||AP|exp(2r 2 /r 2 



(14) 



Figure 5(a) shows the efficiency e for various Rabi fre- 
quencies and pulse delays for a decay rate of T = 0.0 It' 
and a fixed pulse length of T = 800/ The behavior 
described by the adiabaticity requirement is clearly vis- 
ible. For a fixed Rabi frequency the efficiency decreases 
with increasing pulse delay r as predicted by Eq. (fT4| . 
The improvement due to increased Rabi frequency is also 
evident. However, for large Rabi frequencies the pulse vi- 
olates the resolved sideband condition fi\ 7 (i) <C v, leading 
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FIG. 5: (Color online) (a) The efficiency for SCRAP for dif- 
ferent Rabi frequencies f2 and delay times r for the A-system 
in Fig. [2] The pulse width is T = 800/f and the pump and 
Stokes pulse detunings are = AJ = 100^. (b) The effi- 
ciency for a delay r = 320/V. Other parameters are as in Fig. 

Wi 



to a sudden deterioration of the efficiency e at large in- 
tensities, see Fig. 5(b) For small Rabi frequencies the 
efficiency e strongly depends on the delay times. For 
t < 80/V the method is not robust. As the effective de- 
cay rate increases with increasing Rabi frequencies, the 
efficiency slowly degrades for greater laser intensity. This 
leads to a ridge in the t-Q diagram. Because the ef- 
ficiency depends on the pulse delay and the Rabi fre- 
quency, accurate knowledge of these pulse parameters is 
required. This can be moderated by increasing the de- 
tuning of the laser pulses. Because the constraints on 
the pulse length are less severe than for STIRAP the 
population transfer can be faster with SCRAP. Together 
with the far detuning this leads to a improved robustness 
against the detrimental effect of spontaneous decay (49j . 



Stokes pulses (r = 0) with one laser having a fre- 
quency chirp. This is a Raman chirped adiabatic pas- 
sage, sometimes called RCAP (28|, though here we refer 
to it as chirped adiabatic rapid passage (in a A-system) or 
CARP. With this system the Stark shifts are eliminated, 
and the detuning Aj(i) reads 



Aj(t) = 6j(t) = a{t - {J max - J)f). 



(15) 



In order to ensure the adiabatic evolution the chirp rate 
a needs to fulfill \a\ < f2j and |a| T 2 > 1 [with the two- 
photon Rabi frequency ttj given by Eq. (|12|)]. These con- 
ditions arise from a Landau-Zener adiabaticity and from 
requiring completion of a Landau-Zener transfer within 
the time-scale of the pulse. Because there is no limit on 
the pulse duration arising directly from the adiabaticity 
requirements the transition can be fast: it is only limited 
by the narrow bandwidth condition vT ^> 1 [301 ] . This 
in turn reduces the susceptibility to spontaneous emis- 
sion. The resolved sideband condition [4l[ requires that 
v ^> f2j(i), which is easy to satisfy since the system is 
in the far-detuned limit Aj 3> O - Under these condi- 
tions the population transferred to the target state can 
be estimated with the Landau-Zener formula [5(1 [EH 



Pj 



i(oo) = P 7 .„(-oo)(l - e 



-7rAj 



(16) 



where Aj = O. 7 (0) 2 /2|a| = r] 2 n*/(2(A P j) 2 \a\). This be- 
havior is confirmed by our numerical simulation [see Fig. 
6(a)| . It shows the simulated efficiency for different Rabi 
frequencies and chirp rates for a system initially in a state 
with P 0i0 (-oo) = 0.3 and Pi, (-oo) = 0.7. The effi- 
ciency increases rapidly with increasing Rabi frequency 
Do until it reaches a plateau. In this region the transfer 
efficiency is well above 98%. For small Rabi frequencies 
the efficiency e deteriorates with increasing chirp rate. 
However, this effect diminishes for large pump intensities. 
Similarly to SCRAP, the increase in the effective decay 
rate for high Rabi frequencies leads to a slow degrada- 
tion of the transfer efficiency for large pulse intensities. 
For very high laser intensities the efficiency rapidly drops 
due to the violation of the resolved sideband condition 
fij(0) <C v. In a different parameter regime, that of 
small chirp rates a <C 10 -5 ^ 2 , the efficiency oscillates 
with changing Rabi frequency, see Fig. |6(b)| Here the 
evolution is governed by Rabi oscillations between the 
two rotational states which together with the finite pulse 
length leads to large fluctuations in the efficiency. In this 
regime a very precise control of the laser pulses is re- 
quired, so CARP is not robust for very small chirp rates. 
For large Rabi frequencies, the efficiency is essentially in- 
dependent of the chirp rate and CARP provides the best 
robustness against uncertainties in the pulse parameters. 



C. Chirped adiabatic rapid passage 



D. Comparison of the adiabatic passage methods 



Another way of achieving an adiabatic population 
transfer is the application of simultaneous pump and 



Both of the adiabatic rapid passage methods, SCRAP 
and CARP, provide fast population transfer. Together 
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FIG. 7: (Color online) The population transfer for the A- 
configuration of Fig. [2] for different decay rates r = Fj = 
T u and the three methods we consider. The parameters are 
optimised at Tj = r„ = O.Olu for each of the three methods 
and are listed separately in the following. For STIRAP (solid 
line): f2 = u/20, T = 4500 /v, r = 3500/V, and A? = Af = 
0. For SCRAP (dashed line): Q = 7.5u, T = 800/V, r = 
320/V, and A? = Af = lOOv. For CARP (dotted line): fi = 
5u, T = 800/V, r = 0, A? = WOp, At = A? - at with a = 
4.69 x 10~ 5 v 2 . Other parameters which are fixed for all three 
cases are: J max = 1, n max = 2, 77 = 0.1, Po,o(— oo) = 0.3, and 
Pi.o(-oo) = 0.7. 



FIG. 6: (Color online) (a) The efficiency of CARP for different 
Rabi frequencies Q and chirp rates a for the A-system in Fig. 
|U The delay is r = 0, the pulse width is T = 800/^, and the 
Stokes pulse detuning is given by A^ = A^ — at. (b) The 
efficiency as a function of Rabi frequency flo for a small chirp 
rate equal to a = 8 x 10 -6 ^ 2 . Other parameters are as in Fig. 
[3(a)] with pump pulse detuning AJ = 100v. 



with the large detuning of the Raman transition from the 
excited state they provide robustness against the adverse 
effect of spontaneous emission. This is clearly evident in 
Fig. [7] where the efficiency of the three methods is plotted 
against the decay rate. In this figure, the values of the 
parameters Oo, T and r are optimised for V — 0.01 v and 
the detunings given. (In the case of CARP, the value of a 
is also optimal.) The same optimal values have been used 
for fixed parameters in Figs. The parameter V has 
been fixed to a reasonable value for diatomic molecules 
like or CO + . This optimisation gives a fairly wide 
range of parameters in Fig. [7] 

Fig.[JJalso shows that, even though STIRAP is efficient 
for small decay rates, it decreases rapidly for larger de- 
cay rates. SCRAP and CARP are far more efficient and 
their suppression of the excited-state population exceeds 
that of STIRAP within the limits imposed by the exper- 
imental requirements. Even though SCRAP and CARP 
exhibit similar efficiencies CARP is more robust against 
uncertainties in the laser pulse parameters. Hence, we 
use CARP for our investigation of the population trans- 



fer in a multilevel system with J ma x 
Section. 



5 in the next 



IV. CHIRPED ADIABATIC RAPID PASSAGE 
IN A MULTILEVEL SYSTEM (J max > 1) 

Having derived the conditions for efficient population 
transfer with chirped two-photon Raman transitions, we 
turn now to systems where we include a larger number 
of rotational states in the calculation (J ma x > !)■ Con- 
sequently the mapping process involves multiple (J ma x) 
pairs of laser pulses. We investigate the state mapping 
given in Eq. (J5J) , for a thermal initial state distribution 
with (3B = 0.15. This corresponds to a system where 
approximately six rotational levels are significantly pop- 
ulated and consequently we will take J max = 5. In our 
simulation we use a Lamb-Dicke parameter of r\ = 0.1 
and a decay rate of T = 0.0 \v which for realistic trap fre- 
quencies (of the order of a few MHz) corresponds to typi- 
cal decay rates of electronically excited states of diatomic 
molecules. The pump laser is detuned by Aj = 100i^, and 
the Stokes laser is chirped with the rate a = 16 • 10~ 5 ^ 2 . 
Both lasers have a peak Rabi frequency of f2 = 5^, a 
pulse length of T = 800/z/ and a delay between succes- 
sive pulse pairs of f = 4800/za These parameters were 
chosen on the basis of the simulations with two rotational 
levels (Section lllll) and realistic experimental parameters. 



In Fig. 8(a) the initial population distribution over all 
states \J,n) is plotted, along with the final distribution 



(a) Initial population distribution. 




(b) Final population distribution. 

FIG. 8: (Color online) The initial thermal distribution (a) 
for a system with Jmax = 5 and /3B — 0.15, and the fi- 
nal population distribution (b) after the completion of five 
sequential chirped two-photon Raman transitions. The de- 
tuning for the pump pulse is equal: Aj = 100^, the Rabi 
frequency is £l T = hv, and all the decay rates are equal 
to Tj — r„ = O.OIza Other parameters: r\ = 0.1, r = 0, 
f = 4800/V, a = 16 x 10~ 5 ^ 2 , T = 800/V and n max = 6. 



in Fig. 8(b) Apart from some small population loss, 
see Fig. |9J and some weak scattering of population into 
states other than the states |0,n), the two distributions 
agree very well. The total population transferred into the 
rotational ground-state is 92%, while the total population 
loss into the uncoupled states is only 1.5%. 

The remaining 6.5% of the total population is mostly 
left in the initial states \ J > 0,0) [sec Fig. |8(b)"| . Losses 
due to spontaneous emission can be further reduced by 
increasing the detuning. The population which remained 
in the higher, coupled rotational states |J > 0,0) can 
be transferred into the ground-state by sympathetically 
cooling the molecule's motion and reapplying the cool- 
ing pulse sequence. Using this approach the popula- 
tion of the ground-state can be even further increased 
from 92% to above 98.4% with a total loss into the 



FIG. 9: (Color online) The (small) population loss as a func- 
tion of time for the two-photon chirped Raman passage in 
a system with Jm ax = 5, see Fig. [S] The detuning for the 
pump pulse is equal to A^ = lOOf and the Rabi frequency is 
fio = 5f. Other parameters are as in Fig. [S] 



uncoupled state of only 1.6%. Additional simulations 
with Aj = 200v and 10 3 ^, resulted in higher efficien- 
cies, 95.7% and 98.9%, respectively The corresponding 
losses are less than 0.9% for A^ = 200^ and 0.2% for 
AP = 10V 

For large numbers of populated states the transfer ef- 
ficiency can be estimated with the following equation: 



Pi 



al ~ ^ ^ 

,7=0 



P(J-n 



jy 



(17) 



where e is the transfer efficiency for the simple A-system, 
and P(J) is the initial population distribution ([5]). This 
agrees well with our numerical simulation of six rota- 
tional levels. 

Starting from a vibrationally cold system, these val- 
ues of the efficiency e show that we can reach the mo- 
tional and rotational ground-state of molecules [52j . For 
a CO + ion with a rotational temperature of T « 100K, 
the first 15 rotational levels are significantly occupied 
initially with a ground-state population of only 3%. By 
applying the CARP state mapping with a detuning of 
Aj = 100^ and taking the selection rule A J = 0, ±2 
into account the population of the two lowest lying states 
can be can be increased to 85%. Using a detuning of 
Aj = 10 3 ^ this can be improved to 97.8% for a single 
cooling cycle. For a trap frequency of 4 MHz this cooling 
cycle will be completed within 10 ms. 



V. CONCLUSION 

In this work we presented an efficient method to cool 
the internal states of molecules by means of coherent pro- 
cesses (and sympathetic cooling) thus suppressing the 
problematic spontaneous decay into uncoupled states. 
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By coupling the internal molecular state to the motion of 
the molecule, that internal state can be mapped onto a 
motional state. Utilizing this, the internal state is cooled 
close to its ground-state if the molecule's motion was ini- 
tially reduced to the motional ground-state through sym- 
pathetic cooling. Ultimately all the degrees of freedom of 
the molecule can be cooled by the application of sympa- 
thetic cooling to the final motional excitation. Due to its 
high efficiency the method presented here is not only use- 
ful to cool the internal state, but can also be employed to 
detect the internal state of the molecule by measuring its 
motional state with an atom which is trapped alongside. 

We have studied various adiabatic methods for a range 
of laser pulse parameters which are relevant for an ex- 
perimental implementation of this cooling scheme. The 
motion of the ion imposes restrictions on the dynamics 
of the population transfer process which severely limit 
the possible parameter range for the laser pulses. For 
the near-resonant method (STIRAP), population trans- 
fer efficiency is very low accompanied with a large pop- 
ulation of the excited state. Population losses can be 
suppressed, if far-detuned chirped adiabatic two-photon 
Raman passage methods are employed. Schemes that use 
chirped laser pulses (CARP), or self-induced adiabatic 
passage (SCRAP/SIARP) by Stark shifting the transi- 
tion frequencies, turned out to be very efficient. When it 
comes to the comparison of CARP and SCRAP the for- 
mer method has the advantage of easy optimization since 
it has no dependence with respect to pulse shape. Fur- 
thermore, for both methods, and unlike STIRAP, the re- 
solved sideband condition imposes less severe constraints 
on the useful parameter space. 

The requirements for all three methods were derived 
with simulations for a A-system. Using the results from 
this simple model, we were able to demonstrate the ap- 
plicability of CARP in systems with more than two rota- 
tional states. For far-detuned transitions, a high-fidelity 
population mapping from the internal to the motional 
degrees of freedom is possible. Losses were very low and 



our simulations indicate that the fidelity can be further 
improved by detuning the laser pulses further from the 
transition. 

In the scheme we propose here, each rotational level is 
coupled to the excited state by a laser. Due to the large 
rotational level splitting of light molecules this means in 
turn that multiple lasers are required. The number of 
levels N which have a population larger than the cut- 
off population P, and therefore the number of required 
lasers, can be estimated as N « yj— ln(2P)fcsT/S for 
small P. Even though the number of levels N for 
molecules at room temperature can be of the order of 
25, for temperatures of a few Kelvin this reduces to well 
below 10 states. In many experiments molecular ions can 
be prepared in low lying rotational states by employing 
photo-association or state selective photo-ionisation in 
conjunction with supersonic beam expansion or buffer- 
gas cooling. However, due to the interaction with black- 
body radiation and collisions the internal temperature 
quickly thcrmaliscs. By applying the scheme proposed 
here, this thermalisation can be suppressed to main- 
tain the ground-state population. Additionally, the state 
mapping can be employed to detect the internal states of 
the molecule in a non-destructive manner which is bene- 
ficial for high-resolution spectroscopy of molecules. 

In conclusion, we have developed a fast scheme for 
cooling the internal states of single molecules by employ- 
ing adiabatic passage methods which provide a high effi- 
ciency in conjunction with robustness against variations 
in the parameters of the involved laser pulses. 
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